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ABSTRACT 


w 

Interior  ballistic  equations  suitable  for  computer  analysis  are  de¬ 
veloped  to  predict  the  performance  of  the  40mm  grenade  launcher  M79 
firing  the  multiple  projectile  cartridge  XM576.  These  equations  take 
into  account  the  effects  of  the  nozzle -start  and  shot-start  pressures 
which  are  characteristic  of  this  system.  Also  considered  is  the  effect 
of  secondary  igniter  material  upon  the  ballistic  performance.  An  iso¬ 
thermal  model  is  used  for  the  period  during  which  the  propellant  is 
burning.  Only  propellant  of  constant  burning  surface  area  is  considered 
in  this,  study. 


GLOSSARY  OF  SYMBOLS 

A  Bore  Area 

B  Exponential  Burning  Rate  Coefficient 

B'  Linear  Burning  Rate  Coefficient 

C  Propellant  Weight 

D  Propellant  Web 

g  acceleration  Due  to  Gravity 

n  Burning  Rate  Exponent 

Nj  Fraction  of  Charge  Burnt  Which  Remains 
in  High  Pressure  Chamber 

N2  Fraction  of  Charge  Burnt  Which  is  Dis¬ 
charged  Into  Low  Pressure  Chamber 

P}  Pressure  in  High  Pressure  Chamber 

P2  Pressure  in  Low  Pressure  Chamber 

P„_  Nozzle -Start  Pressure 

ns 

P88  Shot -Start  Pressure 

3  Total  Throat  Area  of  The  Nozzles 

t  Time 

tjj  Time  at  Which  Charge  is  Burnt 

tQ  Time  at  Which  Nozzles  First  Start  to 
Discharge 

Uj  High  Pressure  Chamber  Volume 

U2  Low  Pressure  Chamber  Volume 


in.  /sec-(Psi)n 
in.  /  sec -Psi 
lbs 
in. 

ft/sec^ 

Dimensionless 

Dimensionless 

Dimensionless 


psi 


psi 

psi 

psi 

.  2 

in. 

seconds 

seconds 

seconds 


iii 


V . 


GLOSSARY  OF  SYMBOLS  (Coat'd) 


V 

Velocity  of  Projectile 

in.  /  sec 

VB 

Velocity  at  All-Burnt 

in.  /  sec 

W 

Projectile  Weight 

lbs 

X 

Travel 

in. 

*B 

Travel  at  All-Burnt 

in. 

y 

Ratio  of  Specific  Heats 

Dimensionle  8  s 

C 

Heat  Loss  Coefficient 

Dimensionless 

6 

Propellant  Density 

lbs/in.  3 

V 

Co-Volume  of  Propellant 

in.  3/lb 

X 

Effective  Mean  Impetus 

ft-lbs/lb 

*o 

Loading  Density 

gms/cm3 

$ 

y 1/2  [2/<r  +  l)fr  +  lV2,y  ■ 11 

Dimensionless 

* 

t*/Vx7g)(SD/2B'C) 

Dimensionless 

4 

Fraction  of  the  Charge  Burnt 

Dimensionless 

tjto 

Fraction  of  the  Charge  Which  is  Burnt  at 
the  Time  the  Nozzles  First  Start  to 
Discharge 

Dimensionless 

Fraction  of  the  Charge  Which  is  Burnt  at 
the  Time  the  Projectile  First  Starts  to 

Move 

Dimensionless 

*i 

Effective  Fractional  Charge  Increase 
from  Burning  Igniter  Material 

Dimensionless 

iv 


GLOSSARY  OF  SYMBOLS  (Cont'd) 


Ui/6(y  -  1)  S  ^Ag 


Nondimensional  Time  Factor  Used  to 
Describe  the  Period  After  All-Burnt 


seconds 


Dimensionless 
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INTRODUCTION 


The  interior  ballistic  system  for  the  40mm  grenade  launcher  firing 
a  multishot  antipersonnel  round  is  unique.  It  differs  from  that  used  in 
conventional  small  arms  in  that  it  is  a  high-low  ballistic  system  opera- 
ting  with  shot-start  and  nozzle-start  mechanisms.  "Nozzle-start"  means 
that  the  gas  does  not  discharge  from  the  high  pressure  chamber  of  the 
system  until  a  blow-out  disc  which  seals  the  nozzles  ruptures.  This 
does  not  occur  until  a  certain  pressure  (nozzle-start  pressure,  Pn8) 
is  reached  and  a  certain  fraction  (40)  of  the  charge  is  burnt.  Also,  the 
"shot-start"  condition  is  not  reached  until  a  certain  pressure  ("shot- 
start"  pressure  P88)  is  obtained  in  a  low  pressure  chamber.  Likewise, 
this  does  not  occur  until  a  certain  fraction  (%8)  of  the  charge  is  burnt. 

A  shot-start  pressure  is  obtained  by  simply  crimping  the  projectile  in 
the  cartridge  case. 

A  computer  study  has  been  conducted  to  determine  how  changes 
in  nozzle-start  pressure  and  shot-start  pressure  together  with  changes 
in  nozzle  throat  area  and  charge  weight  affect  the  performance  of  the 
system.  This  study  will  appear  in  a  separate  Frankford  Arsenal  re¬ 
port.  ** 

The  interior  ballistics  for  a  high-low  gun  have  already  been  in¬ 
vestigated  by  Corner^  and  Kapur 3,  However,  these  analyses  are  not 
adequate  for  the  system  described  in  this  report  since  they  do  not  take 
into  account  the  nozzle -start  and  shot-start  pressures  and  also  the  ig¬ 
niter  material.  Since  it  was  felt  that  these  factors  would  significantly 
affect  peak  pressures  and  velocity  they  were,  therefore,  taken  into  ac¬ 
count  in  this  study.  Also,  refrences  2  and  3  both  use  a  term  which  rep¬ 
resents  the  linear  burning  rate  constant  for  web  regression.**  This 
quantity  is  most  difficult  to  obtain  without  testing  the  system  experi¬ 
mentally  and  measuring  the  maximum  pressure.  In  this  report  an 
iterative  method  is  presented  for  determining  $  or  B'  from  a  know¬ 
ledge  of  the  gun  and  propellant  parameters. 


See  references. 

jVjV 

In  general  there  are  two  burning  rate  constants  used  in  ballistics. 

($  is  used  to  determine  the  rate  at  which  the  web  burning  from  both 
sides  is  diminished,  and  B1  as  used  in  this  report  to  show  the  rate  at 
which  the  web  burning  from  one  side  is  diminished  {i.e,  B'  •=  ,?/2). 


Numerical  integration  of  the  nonlinear  differential  equations  in¬ 
volved  in  the  interior  ballistics  of  this  system  was  performed  on  a  digi¬ 
tal  comput'r  using  a  Runge-Kutta  method.  * 


ASSUMPTIONS 


1.  The  temperature  of  the  gases  during  the  burning  period  was 
assumed  constant.  This  is  equivalent  to  saying  that  the  impetus  (X)  is 
constant  during  the  entire  burning  period.  This  applies  to  the  gas  in 
both  the  high  and  low  chambers  and  in  the  bore  of  the  gun.  For  numeri¬ 
cal  analysis^  a  value  of  X  equal  to  0.  90  of  the  true  propellant  impetus 
was  used.  This  assumption  was  necessary  to  minimize  the  difficulties 
encountered  in  obtaining  a  solution. 


2.  The  surface  area  of  the  propellant  is  assumed  to  remain  con¬ 
stant  during  the  burning  period.  This  is  equivalent  to  assuming  the  pro¬ 
pellant  form  factor^’  ^  is  zero  and  that  no  unburnt  propellant  is  ejected 
to  the  low  pressure  chamber. 


3.  The  ratio  of  the  pressure  in  the  low  pressure  chamber  (P^) 
to  that  in  the  high  pressure  chamber  Pj  is  assumed  to  be  always  less 
than  |2/(y  +  i )JV /V  -  1  (i,e,  p^/p^  <  |2/(y  +  1)1^^  "  *  ).  Thus  the  flow 
from’the  high  pressure  chamber  is  assumed  to  be  sonic. 

4.  Covolume  effects  on  the  flow  between  the  high  and  low  pres¬ 
sure  chambers  are  neglected, 

5.  The  total  heat  loss  up  to  any  time  t  after  shot-start  is  pro¬ 
portional  to  the  kinetic  energy  of  the  projectile  at  that  time.  C  is  used 
to  denote  the  constant  of  proportionality. 

6.  There  exists  a  uniform  gas  temperature  and  pressure  through¬ 
out  the  volume  behind  the  projectile. 

7.  Flow  in  the  high  pressure  chamber  is  assumed  to  be  quasi¬ 
steady  so  that  no  pressure  waves  or  gradients  are  introduced  into  the 
system. 

8.  Propellant  ignition  is  assumed  isochoric. 
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THEORY 


Derivation  of  Equations 


Internal  Ballistics  of  Gun  Before  All-Burnt 


1.  Equation  of  State  For  High  Pressure  Chamber 


12  CXN, 

P,  =  - - 

1  V!  -  CU  -  4)/6  -  CN jt/ 


2.  Equation  of  State  For  Low  Pressure  Chamber 


P2  = 


12  CXN2 

V2  +  Ax  -  CN2T? 


3.  Equation  of  Web  Regression 


D  =  2B'P, 


4.  Equation  of  Continuity  For  High  Pressure  Chamber 


T  di 


*SP, 


dt  dt  c 


From  the  equation  of  web  regression  (eq  3)  substitute  for  Pj  in 
the  above  equation  (4a) 


dNj  \bSD  d* 

dt  dt  2  B’  c  dt 


I 


the  flow  factor  4  is  defined  as 


4>SD 

2B'  /XTg" C 


Substituting  for  4  in  (4b): 


dNj 

”dT 


(4c) 


(4d) 


Integrating  using  the  boundary  4  =  40  when  =  N0  =  40 


Ni 


=  *  [1  -  *(«  .  «„)/«] 


(4e) 


It  should  be  noted  at  this  point  that  4©  represents  the  fraction  of 
the  charge  which  is  burnt  at  the  time  when  the  nozzle  first  starts  to  dis¬ 
charge.  In  "Appendix  A"  it  is  shown  that  it  is  possible  to  take  into  ac¬ 
count  the  effect  of  the  igniter  material  (igniter  weight  -  Cj ,  impetus  - 
Xi)  upon  the  performance  of  the  gun.  To  do  this  it  is  necessary  that  the 
propellant  weight  (Cp)  be  suitably  adjusted  and  a  determined  based 
upon  die  assumption  that  the  total  igniter  charge  is  burnt  before  any  of 
the  propellant  starts  burning. 

5.  Equation  of  Continuity  For  Low  Pressure  Chamber 

N2  =  4  -  N i  .  4  («  -  4o)  (5) 

6.  Pressure -Time  Curve  in  High  Pressure  Chamber 

Before  All-Burnt 

a.  Derivation  of  Pres  sure -Time  Equation  for  High 

Pressure  Chamber  From  Ignition  to  Nozzle-Start, 

<*  =  *o) 

Assuming  that  4=0  for  the  period  of  ignition  to 
nozzle  start  the  pressure  in  the  high  pressure  chamber  is  given  by 


4 


D  d4 
2B'  dt 


(6a) 


r>  _  12  CA  $ 

1  "  Vx  -  C(1  -*)/6  -  Cr,i 

or 

d$  _  _ 24  B'CX» _ 

dt  D Ijj  -  C/6  +  MC/6  -  CH))‘ 


dt  yt+  T* 

where 

€  =  24  B'CX/D 

T  =  C/6  -  Crj 
jl  *  Vj  ■  C/6 


(6b) 


(6c) 


(6d) 


Integrating  (6c)  with  the  initial  condition  $  =  when  t  =  0 
the  following  is  obtained 


r/e(*  -  $i)  +.i/c  ln(^/4i)  =  t  (6e) 

0MMo 

where  is  determined  by  the  amount  of  igniter  material  burned  initi¬ 
ally  (see  Appendix  A).  Having  thus  determined  $  as  a  function  of  time 
it  is  then  possible  using  (6a)  to  determine  Pi  as  a  function  of  time  for 
this  period. 

b.  Derivation  of  Pres  sure -Time  Equation  for  High 
Pressure  Chamber  From  Nozzle-Start  Until 
All -Burnt 

Substituting  equation  (4e)  for  Ni  in  equation  (1)  the 
following  equation  is  obtained: 


5 


3  -?* 


Pi  « 


12CX^[l  -  4(4-  *0)/*] 


Ui  -  C(1  -  4)/6  -  c»)4[i  -  4(4  -  40)/*J 


substituting  for  Pj  in  equation  (3) 

__D_d*  .  12 CX*[l  -*<4  -»„>/{] 

2B  dt  Uj  -  C(1  -  4)/6  -  Cr*[i  -  4  (4  -  40)/«] 


d4  _  41  +  g*4 
dt  “  +  X'4 

where 

4'  =  24B'CA*04/D 

c'  =  24B'cA(1  -  4 )/  D 
jl'  =  Uj  -  C/6  -  C?j44c 
r'  =  C/6  -  C V  +’  C rj4 


Integrating  (6h)  using  the  initial  conditions  4  a  *0  at  t  =  t0 
the  following  is  obtained: 


where 


where 


r  '/<•  (<1.  _  *0)  +  K/f1  ln( i'  +  <'«/?'  +  f'v.)  =  t  -  t0  (6i) 


‘o  s‘  s  ‘b 


K  =  (f'jl'  - 


& 


< 


and  tb  is  the  time  to  all-burnt  (4  =  1).  Having  thus  determined  ^  as 
a  function  of  time  it  is  then  possible  using  (6f)  to  determine  P}  as  a 
function  of  time  for  this  period  of  time. 

7.  Equation  of  Motion 


12  g  dx 


ap2 


(7a) 


now 


or 


2  B*Pj 


V  = 


2B'P1/dx 


3  1  /  dx  \ 

D  \dw 


dV  _  2  B  d  ,  p 
d5  '  D  d5 


dV 

dV  d4 

dV 

dx 

=  d*  to 

=  d£ 

(7c)  into  (7d) 

dV 

2  B'Pi 

2  B* 

dx 

"  DV 

D 

■ 

■is)] 

p) 


4B^Pi  d 

d* 


d2v 


“  ^  (Pl  5?) 

M 


substituting  (7e)  into  (7a) 


d 

d* 


,  dx 

Id* 


_  3ADzg 
B12W 


© 


(7b) 


(7c) 


(7d) 


(7e) 


7 


or 


d2»  .  3AD2g  /P2\ 

1  d*2  B1ZW  'P1 ' 


(7fl 


since 

P i  =  P,^.*)  aad  Nt  =  Ni(4) 

now 


21  s/8Pn  +  /dpn 

d4  VdNi/tdT  Vat  /N 


from  (4e) 


using  (1) 


likewise 


.  dpl  _ 
•  dT  " 


Nj  =  $[l 

1 

& 

1 

O 

dNi/d*  = 

1  -  * 

/dpl\  _ 

12C\[ui  -  C/6  (1 

-*>] 

\8Ni4 

[U!  -  C/6  (1  -4)  - 

CNilj]2 

(tlL)  = 

12  CX  C/6  {*[l 

-  *  (4  -  *„)/*]  } 

[Ui  -  C/6  (X  . 

.  $)  -  CNj  tj]^ 

5, 

*■< 

o 

CM 

.  C/6  (1  -  4)][l  -  * 

]  -  C«/6[l  -  *(*  -  «0)/«]} 

{Ul  • 

■  C/6  (1  -  4»)  -  CT)<* 

H* 

1 

1 

A 

A 

1 _ 1 

(7g) 

(7h) 


(7i) 


(7j) 


8 


) 

* 


thus  (7f)  reduces  to 


12 CX  {[Ui  -  C/l  -  *)]{1  -  4)  -  C*/6[l  -  ¥(*  -  #0)/«]} 
(Ul  -  C/6  (1  -  *)  -  C»?*[l  -  *  (*  -  *<>)/«]) 


/d2x\  3  AD2g  /  PZ  \ 

\d42y  B12W2  \P1/ 


where  equations  (1)  and  (2)  are  now  substituted  for  Pj  and  P2. 

This  yields  a  nonlinear  differential  equation  which  is  inte¬ 
grated  numerically  using  the  method  of  Runge-Kutta.  *  The  boundary 
conditions  are  dx/d$  =0,  x  =  0  when  4-  =  «8 .  From  this  equation  it  is 
then  possible  to  determine  the  travel  and  velocity  as  functions  of  4- 
and  thence  using  equation  (6i)  as  functions  of  time. 

8.  Equation  for  Maximum  Pressure  in  High  Pressure 
Chamber 

The  primary  difficulty  in  determining  peak  pressure 
from  equation  (7j)  is  that  4  is  a  function  of  B1  which  is  in  itself  a  func¬ 
tion  of  peak  pressure. 


0  SD 


fiu  2B'C 


(8a) 


48  is  the  fraction  of  the  charge  which  is  burnt  at  the  time  the  projectile 
first  starts  to  move.  At  this  time  the  pressure  is  equal  to  the  shot- 
start  pressure  (P88).  4S  may  be  calculated  from  equation  (2)  by  sub¬ 
stituting  P88  —  P2  and  4>s  =  4-  and  letting  x  =  0.  Thus  4g  is  given  by 


Pss  U2 

4  _  d  +  _ 88-  Z _ 

8  ’  0  12  CX  4  +  P88Ct?4 


9 


* 


B1  is  usually  determined  by  using  the  value  of  peak  pres¬ 
sure  obtained  from  experimental  data,  or  the  peak  pressure  is  estimated 
and  then  B(  is  obtained* 

One  way  of  determining'  B*  from  the  peak  pressure  data 
obtained  experimentally  is  to  select  that  value  of  B*  such  that  the  area 
under  the  assumed  linear  burning  rate  (r1  =  B'P)  curve  equals  that  under 
fiifr  exponential  burning  rate  (r  =  BP11)  curve*  The  limits  are  taken  from 
P  s  0  to  P  =  Pmax* 


Thus 


I. 


This  gives 


max  /*rmax 

B'PdP  =  I  BPndP 
o  Jo 


B  = 


2B 


(n  1  1)  P<i«n) 


(8b) 


Values  of  B  and  n  can  usually  be  found  in  literature  or  obtained  from 
burning  rate  studies  for  different  propellants. 

From  equation  (7j)  it  may  be  seen  that  P}  increases 
manotonically  with  4  (i,e. ,  dPi/d#  *  0),  provided  that  the  following 
condition  is  met: 

Ux  (1  -  ¥)  -  C/6  jl  -  *(1  -  *0]  s  0 


With  this  condition  Pi  will  be  maximum  when  $  is  maximum*  This 
will  occur  at  all-burnt  (4  =  1).  Since  loading  density  in  gm/cm^  is: 


60  =  27.7  C/Ui 

the  condition  for  maximum  pressure  to  occur  at  all-burnt  is  thus 
C/U  *  6<l  -  <)/l  -  4(1  -  4-o) 


10 


or 


LQ  *  27.7  C/U  *  27.7  t  (1  -  4)/l  -  4(1  -  * 0 ) 


Some  values  of  4  and  which  appear  feasible  are: 
0.40  <  4  <0.85 
0. 02  <  40  <  0.  20 


Taking  the  case  4  =  0.40  and  40  =  0.20  the  loading 
density  £0  would  then  have  to  be  less  than  about  0.  97  gms/cm^  for 
peak  pressure  to  occur  at  all-burnt.  Thus  our  original  assumption 
that  peak  pressure  occurs  at  all-burnt  appears  to  be  a  reasonable  one 
for  most  practical  high- low  guns  using  propellant  with  a  constant  burn¬ 
ing  surface. 


The  equation  for  the  peak  pressure  is  then  given  by  (1) 

with  4  =  1. 


12  CA  £l  -  4f(l  -  ®0)] 
Ui  -  Cr?fl  -  4(1  -  *o)] 


(8c) 


As  mentioned  previously  4  which  is  dependent  upon  b' 
is  also  a  function  of  Pim.  Thus  an  iterative  technique  was  used  to  de¬ 
termine  Pim»  Pi*8*  a  value  of  4  is  assumed.  This  value  is  then  used 
in  equation  (8c)  to  determine  Pim.  The  value  of  B  is  then  determined 
from  equation  (8b)  and  thence  using  this  value  of  B  a  new  value  of  4 
calculated  from  equation  (8a).  This  is  compared  with  our  original  exti- 
mate  of  4  and  the  procedure  continued  until  both  agree.  The  peak  pres¬ 
sure  is  then  also  established  from  equation  (8c). 


Internal  Ballistics  of  Gun  After  All-Burnt 


For  the  period  after  all-burnt  the  isothermal  model  will  be 
neglected  and  the  gases  in  the  high  pressure  chamber  assumed  to  under¬ 
go  an  adiabatic  expansion.  The  flow  of  the  gas  likewise  will  be  assumed 
to  be  sonic  during  this  period. 


\ 


J 

l 


i 
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9.  Pressure  In  The  -High  Pressure  Chamber  After  All-Burnt 

Following  Corner's  approach  it  may  be  shown**  2  that 
after  all-burnt  the  equation  for  the  pressure  in  the  high  pressure  cham¬ 
ber  as  a  function  of  time  is  given  by: 


12CA[l-4(l-i0j|  -zy/y.l 

Pi  =  - r - T-ll+U-tbVe) 

ux  -T?c[l  -«<1  -*o)J 


Nj:N|  after  all -burnt  and  is  given  by 


nJ  =  [1  -  *<i  -  «„)]  (1  +(t-tb)/e) 


-2/y-i 


g  _ _ » _ 

1  "  12 (y  -  l)*S(Xg)!/z 


10.  To  Determine  the  Muzzle  Velocity  After  All-Burnt 

The  pressure  in  the  low  pressure  chamber  after  all- 
burnt  is  given  by: 


p2  = 


12CXN2  t2 
U2  +  Ax  -  CtjN2 


Since 


1*2  =  (1  -  n;> 


we  have 


p2  = 


12  CX  (1  -  Nj)  T2 
U2  +  Ax  -  CV(1  -  nJ) 


(10a) 
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The  equation  of  motion  has  already  been  given  as: 


W2  V  dV  .  W2  d2x  _  Ap 
12  g  dx  12  g  dt2  2 


(7a) 


The  rate  of  change  in  internal  energy  of  the  gases  in  both  chambers  is 
given  by: 


El  =  -771  £  [NiTl  +  U-w'l)  tz]  <10b> 

where  T j  is  the  ratio  of  the  gas  temperature  in  the  high  pressure  cham¬ 
ber  after  all-burnt  to  the  temperature  which  it  had  during  burning,  and 
likewise  is  the  ratio  of  the  gas  temperature  in  the  low  pressure  cham¬ 
ber  after  all-burnt  to  the  temperature  which  it  had  during  burning. 

The  rate  at  which  work  is  done  on  the  projectile  is  equal 
to  die  rate  of  change  of  kinetic  energy: 


h  “  [w*/dt)J 

24  g  dt  V 


(10c) 


If  it  is  assumed  that  heat  loss  for  any  period  of  time  after 
shot-start  is  proportional  to  kinetic  energy  of  the  projectile  at  that  in¬ 
stant  then  the  rate  of  loss  01  energy  by  heat  transfer  can  be  written  as: 


E3 


C  W  d_ 
24  g  dt 


£(dx/dt)2 


(10d) 


where  C  is  the  proportionality  constant. 

The  sum  of  (10b),  (10c)  and  (lOd)  must  equal  zero  since 
the  total  energy  of  the  system  is  constant. 
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Let  (y  -  1)  =  (  y  .  1)(1  +  C)  in  the  above  sum.  Integrating 
this  stun  and  making  use  of  the  following  initial  conditions  =  N ^ , 
rl  s  Tlb  =  r2  ~  r2b  =  *  subscript  b  indicates  the  condition  at 
all-burnt)  when  (dx/dt)  =  (dx/dt)b  we  obtain: 

[n;  TJ  -  Nlb  +  (1  -  n’,)  r2  -  (1  -  Nlb>]  = 


-  —  jj  dx/dt)2  -  (dx/dt)2  J  (10e) 

Using  equation  (7a),  the  equation  of  state  for  the  low  pres¬ 
sure  chamber  after  all-burnt  (equation  (2))  may  be  rewritten  as: 

12CX(1  -  Nj)  t2  *  P2[o2  +Ax  -  C(I  -  N j ) tj ] 

+  <10£> 


Using  equation  (9b),  the  equation  of  state  for  the  high 
pressure  chamber  after  all-burnt  (equation  (1))  may  be  rewritten  as: 

12 CXn'j  Tj  x  pJuj  .  CNjr)] 

»j]  [i  +  (t  -  v/e]  ylr  (iog) 

substituting  lOf  and  lOg  into  lOe 

pib [u i  -  cm;,]  [i  +  (t-tb)/e]"Zy/y"1  +  0  [u2 

+  Ax  -  C(1  -  N \)ri~\  -  1 2  C X  -  -  [(dx/dt)2 

-  (dx/dt)^] 
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=  plb  [ux  -  cn; 


1 


or  multiplying  through  by  A  and  rearranging  terms: 


■j—  (d2x/dt2)  [U2  +  Ax  -  C ri(l  -  Nj)]  +  (y  -  1 ) AW / 24g  (dx/dt)2 

=  (y  -  l)WA/24g  (dx/dt  +  12ACX  -  Plb[u! 

-2y/y-l 


-  cnJt?]  [i+(t-t^)/e 


(10h) 


where 


Ni  ■  N1B[l+tt-tB)/8]’2/y'1  =  [l  -  4(1  -0o)]|_l+(t-tb)/8j 


n-2/y-l 


and  where 


6  = 


2  Ui 


12 (y  -  1)*S  yrr 


Equation  (lOh)  is  a  nonlinear  differential  equation  which 
is  integrated  numerically*  using  the  method  of  Runge-Kutta.  The 
boundary  conditions  are  x  =  xB  and  dx/dt  =  (dx/dt)^  when  t  =  t^ . 


CONCLUSIONS  AND  RECOMMENDATIONS 


1.  The  assumptions  made  and  the  equations  developed  in  this 
report  appear  reasonable  for  predicting  the  performance  of  a  high- low 
pressure  gun  using  nozzle-start,  shot-start  and  secondary  igniter  ma¬ 
terial. 


2.  Future  studies  should  be  conducted  to  determine  how  the 
equations  developed  in  this  report  could  be  modified  to  further  describe 
high-low  pressure  guns  where  factors  such  as  propellant  loss  between 
the  chambers,  progressiveness  and  degressiveness  of  the  propellant 
end  changes  of  gas  temperature  during  burning  are  important. 


< 


APPENDIX  A 


To  determine  the  effect  of  igniter  material  on  the  Interior  Ballis¬ 
tic  Performance. 

Cp  a  Wt.  of  propellant  (lbs) 

Cj  a  Wt.  of  igniter  (lbs) 

C2  a  Effective  charge  wt  of  igniter  (lbs) 

X  a  Impetus  of  propellant  (ft-lbs/lb) 

a  Impetus  of  igniter  (ft-lbs/lb) 

C  a  Total  effective  charge  wt.  (lbs) 

a  Covolume  of  igniter  material  (in.  /lb) 

Assume  that  the  total  igniter  charge  is  burnt  at  the  time  the  pro¬ 
pellant  first  starts  to  burn.  Let  be  the  fraction  of  the  total  effective 
charge  (C)  burnt  at  this  time. 

^i  =  ^/(Cp  +  C2)  (A-l) 

Let  P^  be  the  pressure  generated  by  the  igniter. 

Pi  =  UCiXi/Ui  -  Cp/6  -  r\  •  q  (A-2) 


Assume  this  igniter  pressure  to  be  generated  by  an  initial  charge 
burnt  with  an  impetus  (X )  equal  to  that  of  the  propellant  alone: 


pi  = 


12  (Cp  +  C2)^i  X 


Ui 


(Cp  +  C2)  (1  -  q) 


-  (Cp  +  c2)  V*i 


(A- 3) 
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Substituting  for  4^  from  equation  (A-l) 


12C2X 

p*  *  W,-Cp/«-  c2» 

Solving  for  C2 

Uj  -  Cp/6 
°2  =  UX/Pj  +  T? 

Thus  using  (A-2)  and  (A- 5)  we  can  then  determine 

*i  *  C2/(C2  +  Cp) 
and  the  total  effective  charge 
C  a  C*2  +  Cp 

Likewise ,  the  total  effective  web  is 


(A-4) 


(A-5) 
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